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The dynamics of a bi- dimensional dense granular packing under cyclic shear is experimentally 
investigated close to the jamming transition. Measurement of multi-point correlation functions are 
produced. The self-intermediate scattering function, displaying slower than exponential relaxation, 
suggests dynamic heterogeneity. Further analysis of four point correlation functions reveal that 
the grain relaxations are strongly correlated and spatially heterogeneous, especially at the time 
scale of the collective rearrangements. Finally, a dynamical correlation length is extracted from 
spatio-temporal pattern of mobility. Our experimental results open the way to a systematic study 
of dynamic correlation functions in granular materials. 

PACS numbers: 64.70.Pf, 05.40.Ca, 45.70Cc, 61.43.Fs 



The dynamical behavior of granular media close to the 
'jamming transition' is very similar to the one of liquids 
close to the glass transition [l[. Indeed granular media 
close to jamminffdisplay a similar dramatic slowing down 
of the dynamics @, Lj as well as other glassy features such 
as aging and memory effect 4] . Recently, a "microscopic" 
confirmation of the above similarity has been obtained 
analyzing directly the grain dynamics under cyclic shear 
during compaction 0] or at constant density 6]. The 
typical trajectories of grains display the so-called cage 
effect and are remarkably similar to the ones observed 
in experiments on colloidal suspension Q and in molec- 
ular dynamics simulations of glass- formers Q. As for 
glass- formers, and contrary to standard critical slowing 
down, this slow glassy dynamics does not seem related to 
a growing static local order. For glass- formers it has been 
shown numerically 0, 0, 0, and experimentally ^| 
that instead the dynamics becomes strongly heteroge- 
neous and dynamic correlations build up when approach- 
ing the glass transition. The existence of a growing dy- 
namic correlation length is very important to reveal some 
kind of criticality associated to the glass transition [lj • 
Here we unveil that also granular materials are strongly 
dynamically correlated close to the jamming transition. 
First, we shall focus on two point functions, in par- 
ticular the self-intermediate scattering function, whose 
slower than exponential relaxation suggests dynamic het- 
erogeneity. Then, following recent theoretical sugges- 
tions llfl^ we shall turn to four point correlation 
functions. They have been introduced for glass- formers 
to measure properly dynamic correlations and in- 
deed reveal that the dynamics is strongly correlated and 
heterogeneous. Finally, we shall focus on spatio-temporal 
pattern of mobility, out of which we extract a direct mea- 
surement of a dynamical lengthscale. Our experimental 
results, to our knowledge the first direct measurement 



of four point spatio-temporal correlation functions [l^ . 
open the way to a systematic study of dynamic corre- 
lation functions in granular material as a way towards a 
better understanding of glassy and "jammy" materials in 
general. 

The experimental setup, a more robust and better de- 
signed version of the one presented in p is as follows : a 
bi-dimensional, bi-disperse granular material, composed 
of about 8.000 metallic cylinders of diameter 5 and 6 
mm in equal proportions, is sheared quasi- statically in 
an horizontal deformable parallelogram. The shear is 
periodic, with an amplitude 9 max = ±5°. The volume 
fraction (<1> = 0.84) is maintained constant by imposing 
the height of the parallelogram. We follow 2818 grains 
located in the center of the device to avoid boundary ef- 
fects with a High Resolution Digital Camera which takes 
a picture each time the system is back to its initial po- 
sition = 0. The unit of time is one cycle, a whole 
experiment lasting 10.000 cycles. The unit of length is 
chosen to be the mean particle diameter d. These con- 
ditions are very similar to 6] and by repeating the same 
analysis we find a cage radius of 0.2 and a cage lifetime 
of 300. As discussed in the diffusion is isotropic, at 
least far from the borders of the experimental cell. 
Let us first focus on the self intermediate scattering func- 
tion which measures the dynamics of single particles: 

F s = (F s (k,t)) = 1 £<exp i-ik( rj (t) - rj (0))]>, (1) 

3 

where rj(t) is the position of the jth particles at time 
t. F s (k,t) denotes the non averaged instantaneous ob- 
servable. (•) means (here and in the following) a time 
average over 300 steps of 10 cycle each computed after 
a few thousand cycles, when the systems has reached a 
steady state (at least on the timescale of the experiment). 
The sum in (pp) is over all tracked particles. 
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FIG. 1: On the left: F s (k,t) as a function of time for differ- 
ent odd values of the wave- vector k = 1, 3, 29 from top to 
bottom (as indicated by the arrow and the increasing k /*). 
The black lines are fits of the form exp[-(t/r(/c)) /3(fc) ]. On 
the right: r(k) (top) and (3(k) (bottom) as a function of k. 

The function F 8 (k,t) is plotted on the left of Fig. [I] 
as a function of time for different odd values of k rang- 
ing from 1 to 29. Contrary to glass- formers there is no 
visible plateau in this correlation function although from 
trajectories it is possible to identify a clear cage effect 
(see fig. 2 of la). Note that the short-time dynamics is 
subdiffusive p probably because granular media don't 
experience thermal relaxation. Therefore the separation 
of timescales and the corresponding plateau in correlation 
functions is much less pronounced as in colloids 7]. An- 
alyzing the curves in Fig. we find that the decreasing of 
F s (k,t) is slower than exponential in time. A good fit is 
provided by a stretched exponential: exp[— (t/r(k))^^]. 
We plot on the right of Fig.Q]r(/c) (top) and f3(k) (bot- 
tom) as a function of k. At small k the relaxation time 
scales as k~ 2 and the exponent f3{k) is one. As expected, 
the grains perform a Brownian motion on large length 
and time scales and therefore F s (k,t) exjp(—Dk 2 t) for 
small k and large t. Increasing k the stretched exponent 
decreases and is of the order of 0.7 for k of the order of 
27r, corresponding to the inter-grain distance, and even 
lower for higher values of k. A very similar behavior has 
been found for glass- formers |3, |l3[ . Also the decrease of 
r(k) steepens and decreases sharply for large k. This is 
also related to the short time sub-diffusive dynamics. In 
this regime the particle displacement distribution has a 
variance scaling as t 1 / 2 (not t like for standard diffusion) 
and is well fitted by a Gaussian. Assuming this func- 
tional form and Fourier transforming to get the interme- 
diate scattering function, one finds that the relaxation 
time goes as & -4 , hence a crossover from a k~ 2 behavior 
at small Ho a more rapid decreasing at large k as in 
Fig.ffl 

Dynamical heterogeneity is one of the possible explana- 
tion of the non-exponential relaxation of F s (k, t): the re- 
laxation becomes slower than exponential because there 
is a strong spatial distribution of timescales 0. This 
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FIG. 2: Q a (t) as a function of time for a = 0.05, 0.1, 0.5. 

effect is strong for intermediate and large values of k. 
Instead for smaller fcs, the heterogeneities are averaged 
out and a larger value of /3(fc), going to one for k — > 0, 
is expected. This coincides indeed with the trend found 
in Fig. U However this is not the only possible scenario 
[ill EH • I n the following we want to go one step further 
and show direct "smoking gun" evidences of dynamical 
correlations. For this purpose it is of interest to consider 
the structural relaxation and not only the single parti- 
cle one, as given by F s (k, t). We focus on the the density 
overlap E3 following previous works on glass-forming liq- 
uids: 

Q a = (Q a ) = ^J drdr'(5p(r,t)w a (r-r')5p(r',0)), (2) 

where p(r,t) = ^^(r - n(t)) and 6p(r,t) = p(r,t) - 
(p). The overlap function is a non- normalized Gaussian: 
iy a (r) = exp(— r 2 /2a 2 ). The evolution of Q a (t) is a mea- 
sure of how long it takes to the systems to decorrelate 
from its density profile at time t = 0. Fig. shows that 
the behavior of Q a (t) is similar to the one of F s (k,i), 
as for glass-formers The proper way to unveil 

spatio-temporal correlations is through the fluctuations 
of the temporal relaxation Those are characterized 
by dynamical susceptibilities: 

X f(M) = N{(H(k,t)-(H(k,t))) 2 ) (3) 

where H can either be F s (k,t) or Q a (t). They unveil 
dynamic correlations exactly as fluctuations of the mag- 
netization unveil magnetic correlations close to a fer- 
romagnetic transition, see e.g. [13, One way to 
understand how such susceptibilities relate to spatial 
heterogeneities of the dynamics is to decompose, say, 
Q a (t) in local contributions: NQ a (t) = pj drq a (r,t) 
where q a (r,t) = ~ f dr f 5p(r,t)w a (r — r f )5p(r f , 0). Using 

this expression one finds x ? W = 9 J drG® (r, t) where 
G$(r,t) = ([q a (r,t) - (q a (r,t))}[q a (0,t) - (q a (0,t))}) is 
the spatial correlation of the local temporal relaxation : 
if at point an event has occurred that leads to a decor- 
relation of the local density over the time scale t, (r, t) 
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FIG. 3: Top:xf s (£) as a function of time for odd val- 
ues of k = 1,3, ...,29. Inset: Log-Log plot for k = 
7,9,11,13. Bottom:xJ(t) as a function of time for val- 
ues of a = 0.05, 0.1, 0.5. Inset: Log-Log plot for a = 
0.1,0.15,0.2,0.25. 

is the probability that a similar event has occurred a dis- 
tance r away, within the same time interval t. xj(£), its 
volume integral, quantifies how much spatially correlated 
is the dynamics (see || 0, |2(1 for a more detailed 
discussion). 

Top of Fig. 01 displays xf s 0) for k = L3, ^,29. It has 
the form found for glass- formers 0, UJl LLJ : ^ * s 
of the order of one at small and large time and displays 
a peak at a time somewhat larger than the timescale of 
the sub-diffusive regime. The peak is a clear signature 
of dynamic heterogeneity and shows that the dynamics 
is maximally correlated on timescales of the order of the 
relaxation time. A rough estimation of the correspond- 
ing dynamical correlation length is obtained identifying 
the peak of xf s W> °f the OI> der of 100, to a correlated 
area 7r£^ et , leading to a length £h et ex 6 in agreement 
with a previous estimate p. We find very similar re- 
sults for x±{t) as shown in the bottom of Fig. [3| for 
a = 0.05, 0.1, 0.5. The largest X±(t) is obtained for 
a = 0.15, which corresponds to the typical displacement 
during the sub-diffusive regime. Large and small values 
of a corresponds to small values of the peak because the 
dynamics on small lengthscales is certainly not very cor- 
related and on very long lengthscales the heterogeneous 
character of the dynamics is averaged out. As discussed 
in [2l|, the power law growth of X4* W ( res P- X?W) De_ 



fore the peak with exponents between 1 and 2/3 (see 
insets of Fig. 3) suggest that the dynamic correlations 
cannot be induced by independent defect or free volume 
diffusion. 

We now focus on spatio-temporal pattern of mobility. 
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FIG. 4: Grey-scale plot of q a s {r,t) : at t = 42,435, 1113, 2526 
from top to bottom (a = 0.15). Black regions correspond to 
lower values of . The displacements of the particles during 
the interval of time t are plotted in yellow. The yellow dots 
are particles that have been lost during tracking. 
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Fig. El presents a grey-scale plot of q%(r,t) = J2i — 
ri(0))w a (ri(t) - n(0)) for t = 42,435,1113,2526 and 
a = 0.15, where S(r) is approximated by a Gaussian 
of width 0.3. By definition q^{r^t) measures a coarse 
grained mobility: if the particle that was close to r at 
t = moved away more than a in the time interval t 
then g"(r, t) ~ 0. The yellow lines in Fig^are the par- 
ticle displacements in the time interval t. At short- times 
(t = 42) only few particles have moved and from FigEI 
it appears that they do so in a string-like fashion. On 
larger times (t = 435, 1113) the relaxed regions are rami- 
fied and finally, at very long time (t = 2526) the majority 
of the particles has moved substantially but there remain 
few (rather large) regions not yet relaxed. These findings 
similar to the ones found in simulation of supercooled liq- 
uids 0, 0, 0, suggest that the mobility is organized 
in clusters, which are the direct visual evidence of the dy- 
namical heterogeneities. Left of Fig. displays the prob- 




0.2 0.4 0.6 0.8 1 



FIG. 5: Left:self part of the Van Hove correlation function 
after angular integration at t = 438; the continuous line is 
the pdf for a Gaussian distribution. Right: ln((r4(r, 438)) as 
a function of r (for a — 0.15); the straight line is a linear fit. 
ability distribution of the grains displacement for t = 438 
(corresponding to the maximum of x±(t)) and quantifies 
the excess of fast and slow grains compared to the Gaus- 
sian distribution (in continuous line). On the right of 
Fig. G4(r, 438), the radial autocorrelation of av- 
eraged over ten realizations (for a = 0.15), exhibits an 
exponential decay over a characteristic dynamical length 
£ = 7, in agreement with the value obtained from the 
peak of Xa s • For comparison, the dynamic lengthscales 
reported in experiments close to the glass transition are 
of the order of 5 — 10 molecular diameters |l3j . 
In conclusion our results, to our knowledge the first ex- 
tensive experimental study of two and four point spatio- 
temporal dynamic correlation functions, furnish a direct 
experimental evidence that granular materials close to 
jamming have an heterogeneous and correlated dynam- 
ics. It would be certainly worth studying the possible re- 
lations between the dynamic correlations we have found 
and the diverging length scales that have been proposed 
to show up at the jamming transition (coming from the 
jammed phase) |2^. Our results reveal a remarkable sim- 
ilarity with glass-forming liquids that reinforces the con- 
nection between glasses and jamming systems 1]. They 
open the way to further analysis, varying a control pa- 



rameter (as is the temperature for liquids and packing 
fraction for colloids), or during compaction. That would 
give other important information on the microscopic dy- 
namics and provide stringent constraints to the theory of 
glassy and "jammy" materials in general. 
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